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Abstract 

The target mass corrections are calculated for all structure functions of neutral and charged 
current deep inelastic scattering in lowest order in the coupling constant. Representations 
of the correction to the twist-2 and twist-3 contributions are derived both in Mellin-n 
and x-space. The impact of the target mass corrections on the general relations between 
the twist-2 and twist-3 parts of the structure functions is studied and three new relations 
between the twist-3 contributions are derived. 
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1 Introduction 



Deep inelastic lepton-nucleon scattering provides one of the cleanest possibilities to study the 
nucleon structure at short distances. The light-cone expansion |]T| proved to be one of the most 
powerful techniques in the Bjorken limit, Q 2 ,P.q — > oo,x = Q 2 /2P.q = const., to derive the 
structure of the scattering cross sections, relations between the structure functions and their 
scaling violations. 

The early unpolarized deep inelastic scattering experiments and most of the experiments, 
in which polarized lepton scattering off a polarized target 0, f| has been studied so far, operated 
in the range of lower values of Q 2 . In this domain nucleon-mass corrections cannot be neglected. 
The target mass effects of 0((M 2 /Q 2 ) k ) form one contribution to the power corrections. Unlike 
the case for dynamical higher twist effects [5-7] the target mass corrections can be calculated in 
closed form in all orders in M 2 /Q 2 for deep inelastic structure functions. Because the product 
of a mass factor M k and a genuine twist-/ operator forms an operator of twist r = k + I the 
individual terms in the Taylor-expansion in M 2 /Q 2 of the target mass corrections mix with the 
dynamical higher twist operators under renormalization. Therefore a fully consistent description 
would require a common treatment of both contributions. Little is known so far on the relative 
strength of dynamical higher twist operators and their scaling violations. Before one may try 
to pursue a common treatment various conceptional problems concerning higher twist operators 
have first to be solved. Due to this we limit the present analysis to a systematic study of the target 
mass corrections for all deep inelastic structure functions extending earlier investigations [§, |J . 

Two methods were proposed in the literature for the evaluation of the target mass corrections. 
Nachtmann [10] translated the usual power-series expansion into an expansion of operators of 
definite spin. The problem may be solved by applying the representation-theory of the Lorentz 



group. A second method is due to Georgi and Politzer 0, see also [PH , in which the individual 
nucleon mass terms are collected and resummed. In the domain of large values of x, x ~ 1, these 
methods lead to different results for the operators of lowest twist. This is the region, however, in 
which also the dynamical higher twist terms of any order contribute, which have to be included 
into the analysis in this kinematic domain as well. 

We will apply the method of Ref. || to all polarized structure functions.^ The contributions 
due to the twist-2 and twist-3 operators are studied individually. Both the corrections for 
the Mellin moments of the different structure functions and the inverse Mellin transform to 
Bjorken-x space, which can be derived analytically, are provided. As will be shown, the complete 
resummation of the different series in (M 2 /Q 2 ) k is required to avoid artificially large terms in 
the region x — ► 1. 

A second goal of our investigation is to study the effect of the target mass corrections on 
the integral relations between the different polarized structure functions in lowest order in the 
coupling constant. As was shown in a previous analysis [ITH] the twist-2 contributions to the 



polarized structure functions are connected by three (integral) relations, the Dicus-relation |L6 



the Wandzura-Wilczek relation IJT7J and a new relation given in Ref. |Tj|. A general relation 
which holds for all spins n could also be derived for the valence part of the twist-3 contributions 
to the structure functions g 2 and g 3 . As in Ref. [T5|] nucleon mass effects were disregarded, 
except of those implied by the kinematics in the Born cross sections, twist-3 contributions to 
the structure functions g±,g4 and g$ were not obtained. This picture has to be regarded as 
partly incomplete, since nucleon mass effects have either to be accounted for thoroughly or to 

2 The target mass corrections for <?i and 172 were derived in |l2|, [l3| and to the polarized Bjorken sum rule 
in El following Ref. pi. 
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be neglected at all. In the latter case, however, the scattering cross sections for longitudinally 
polarized nucleons would even not contain the structure functions g^ and and therefore not 
allow to obtain two of the twist-2 relations, namely the Wandzura-Wilczek relation and the new 
twist-2 relation of Ref. [|15[| . On the other hand, a rather asymmetric picture is yet obtained 
for relations in the twist-3 sector. Particularly, there is no relation yet for the complete twist- 
3 term contained in g% and a potential other structure function. As will be shown below the 
complete inclusion of the target mass corrections solves this problem since all polarized structure 
functions obtain, besides the twist-2 contributions, by virtue of these corrections also twist-3 
terms. Similar to the case of the twist-2 contributions observed before three new relations are 
implied for the twist-3 parts of the structure functions in lowest order in the coupling constant. 
These relations provide the possibility of a thorough test of the twist-3 structure of the nucleon. 

The paper is organized as follows. In section 2 the cross sections are summarized for neutral 
and charged current deep inelastic scattering allowing also for contributions due to current 
non-conservation. In section 3 the structure of the forward Compton amplitude is derived and 
section 4 provides the general expressions of the operator product expansion in the presence of 
target mass corrections. The nucleon matrix elements are evaluated in section 5. The relation 
of the moments of the twist-2 and twist-3 contributions to the different polarized structure 
functions are given in section 6. In section 7 the x-space representations for the twist-2 and 
twist-3 contributions of the target mass corrections to the polarized structure functions are 
provided. Section 8 deals with the effect of the target mass corrections on the twist-2 relations 
between the polarized structure functions in the massless limit. Three new relations are derived 
between the twist-3 contributions of the polarized structure functions in section 9. There also 
other twist-3 relations are discussed. Section 10 contains the conclusions and an appendix deals 
with the quark-mass corrections in the case of the Wandzura-Wilczek relation. 



2 The Scattering Cross Section 



The differential Born cross section for polarized lepton-polarized nucleon scattering is given by 



d 3 a 



,2 



//n Evi(Q 2 )Lrwr . (i) 



dxdydcj) Q 

Here the index % denotes the different current combinations, i.e. % = I7I 2 , \jZ\, \Z\ 2 for neutral 
and i = {W^l 2 for charged current interactions. <fi is the azimuthal angle of the final-state lepton, 
x = Q 2 /(2P.q) = Q 2 /(2z/) and y = P.q/k.P are the Bjorken variables, where q = k — k' the four 
momentum transfer to the hadronic vertex and Q 2 = —q 2 . P, k and k' are the proton-, initial- 
and final state lepton four-momenta, respectively. The factors rji(Q 2 ) denote the ratios of the 
corresponding propagator terms to the photon propagator squared, 

v hl \Q 2 ) = l, 

r} hz\( 2 ) _ G F M 2 Q 2 
V W) ~~ 2V2na Q 2 + M 2 > 

v lzl \Q 2 ) = (v hzl ) 2 (Q 2 ), 

V {Q) ~ 2 \ Ana Q 2 + M 2 W ) ' {Z) 

a is the fine structure constant, Gp the Fermi constant and Mz and My/ are the Z and W boson 
masses. 
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The leptonic tensor has the following form 



= E [^(^'' ^OTmUv + 9Al5)u(k, A)] t A')7„(^ + 9Al5)u(k : A) 



(3) 



Here, A denotes the helicity of the initial state lepton. The indices %\ and %2 refer to the currents 
forming the combinations % in Eq. ([!]), and the vector and axial vector couplings read 



9l 
9 Z v 
9T 



9a 



- 2Q I sm 2 6 w , 9 Z a 



n w ~ 
9 A 



o, 
1 

~2 ' 
-1 . 



(4) 



Here Qi denotes the charge of the initial-state lepton. $w is the weak mixing angle. 
The hadronic tensor is given by 



W, 



1 

Air 



d A xe iqx (PS | [Jj(x) f ,4 2 (0)] | PS) . 



(5) 



S denotes the four- vector of the nucleon spin with S • P — 0. In the following we normalize 
S 2 = — M 2 , where M is the nucleon mass. In the framework of the quark-parton model the 
currents are given by 

= £ 'I'A^M 1 + 9 J A f ^)qf(x)U ff , (6) 
/,/' 

where g^ A are the electroweak couplings of the quark labeled by /. For charged current interac- 
tions Uffi denotes the Cabibbo-Kobayashi-Maskawa matrix and gy = 1, — — 1, whereas for 
neutral current interactions Uff> = Sff, g v = e q ,gA = for 7, and 



9l 



9 q v 



1 4 • 2 a 

- - - sin 9 W , g\ 

1 9 
~2 3 ^ 



2' 
1 



for q = u,c,t 
for q = d,s,b 



(7) 



for Z-boson exchange. 

The hadronic tensor is constructed requiring Lorentz and time-reversal invariance. The 
general structure of the hadronic tensor is 



W, 



fJbV 



q 2 



l2\ , PpPv 



q \pa 

2P-q 



Fi(x, (?) + -f^F 2 (x } Q 2 ) - ie^ 



F 3 (x,Q 2 



+ ^F i ( X ,Q 2 ) + ^p^F 5 ( X ,Q 2 ] 



P-q 



2P-q 



q x S rr 



+ 



P-q 

PflSy + S^P^ 

2 



S ■ q 



P P 
P ■ q 



q x (P -qS a -S- qP a ) 

{PIT 2 

9z{x,Q 2 ) 
P-q 



9i{x,Q 2 



s ' lif^.'/ ;(•'•• Q 2 ) + (-9r + ^) { 4^9^, Q 2 ^ 



(P ■ q) 2 
PaS\ 
P ■ q 
(P^Qu + q^Pu)S ■ q 



q 2 P-q 

+i£^\*y^ge(x, Q 2 ) + S ■ q p 9 " g 7 (x, Q 2 ) 



2(p- q y 



9s{x,Q 2 



S^q v ~\~ S u q^ 
2P-q 



99{x,Q 2 



where 



P-q 

~J 2 



Sit — s,. 



S-q 

—2-9/1 
q l 



Here the current indices were suppressed. In general the hadronic tensor is determined by 
five unpolarized structure functions Fi and nine polarized structure functions (ft. Due to the 
structure of the electroweak couplings g\ A given above and due to current conservation, for 
photon exchange only the structure functions and g\p contribute. For the weak currents in 
general all structure functions given above are present. The notation for the structure functions 
Fi, gi and #2 is widely unique in the literature. However, different notations are used for the 



structure functions (ft^ =3 - We follow the definitions of Ref. |L5| for the structure functions (ft| i=3 - 



For the current non-conserving structure functions our definitions coincide with those introduced 
in Ref. |I§, cf. also 0. 

From Eqs. ([l]|§) and (§) one obtains the differential scattering cross sections of a lepton 
with helicity A off a polarized nucleon. For convenience we will consider two projections of the 
nucleon spin vector, choosing the spin direction longitudinally and transversely to the nucleon 
momentum. In the nucleon rest frame one has 



S L = (0,0,0,M) , 

St = M(0, cos a, sin a, 0) . 



(9) 



The polarized part of the scattering cross section for longitudinal nucleon polarization, integrated 
over the azimuthal angle 0, reads 



d?a(\,±S L ) 
dxdy 



^s—Y^CMQ 2 ) 

, / 2xyM 2 \ , vx 2 M 2 . AxM 2 ( xyM 2 \ , 

-2Xy \2-y- ) xg\ + 8A^— p 2 + — — 1 - y - g\ 



2 1 + 



2xM 2 
S 



S 



i-y 



s 



xyM 2 
S 



g\ - 2xy 2 1 + 



S \ 

2xM 2 



S 



S 



9s 



xyM 2 . / xyM 2 \ ,' 

Correspondingly, for transversely polarized nucleons one obtains 



(10) 



d 3 a{\, ±S T ) 
dxdydcj) 



Q 4 



x 2i 



<M 2 



S \ 



xy 



xyM 2 
S 



cos a 



1 / 2xyM 2 \ • 2 / 

y{ 2 - y -—) 3l+ y{ l - V 



-2\yxg\ - AXxgl 
xyM 2X 



:ir 



S 



9l + 1xyg\ - 2Xg l 6 - g\ 



Here = 1, C^ z = g v + Xg A , C z = (g v + Xg A ) 2 , and C 



;i±a). 



As well-known, the contributions of the structure functions g\ and g\ are suppressed by a 
factor of M 2 / S in the longitudinal spin asymmetries 



A L = d 2 a(\, S L ) - d 2 a(\, -S L ) . 



(12) 
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However, they contribute at the same strength as the other structure functions to the transverse 
spin asymmetries 

A T = d 2 a(X, S T ) - d 2 a{\, -S T ) , (13) 

which on their own behave oc Mj \f~S~. The current non-conserving structure functions g\ and g\ 
do not contribute to the cross sections when the lepton masses are neglected. 



3 The Forward Compton Amplitude 

The forward Compton amplitude T„ v is related to the hadronic tensor by 



W { = — Im T 



where 



T; u = iJ d*xe^(PS\(TJ^(x)J?(0)\PS) . 



(14) 



(15) 



The Compton amplitude can be represented in terms of the amplitudes T\. and A\, which are 
related to the unpolarized and polarized structure functions, as 



rpl 



~Qiiv "1 



q ^ Iftf, u) + ^n(x, Q 2 ) - ie^^Tttf, v) 



+ ^(^-) + ^ M ^?^^(g 2 ,,) 



<l X S° *i,.2 ., , ... q\P-qS° -S-qP", , 



+ie fluX a 1 j^A\(q 2 , v) + ie^ Xa 

P Q _|_ Q p p p 

+ — ~ ^ ^ - S -q-£ — 

P • q 



M 4 



A\{q\v) 



4(gV) 
M 2 



<"" "<- 2 - ' < W ^ {S ' q) Ai(q 2 ,v), 



+S-q^Al(q\v) + l-g, u+ f 



M 2 



+^A^^(g 2 , v) + S. q qj ^A^{q\ *) 
, [P^qv + q^Pu)S ■ q j , 2 , + 5^ • 2 \ 



2M 4 



2M 2 



(16) 



The structure functions Fi(x,Q 2 ) and gi(x,Q 2 ) are obtained from the amplitudes Tj(g 2 ,z/) and 
^i(gV) by 



F!(x,Q 2 ) 

^2,3,4,5^, Q 2 ) = 
for the unpolarized structure functions, and 



— Im T^q 2 ,^) , 



1 v 



Im T 2i3A5 (q 2 ,u) , 



(17) 



^1,3,5,6,9(3:, Q 2 ) = 7T77^ Im ^l,3,5,6,9(gV) » 



#2,4,7,8 (X, Q 2 ) 



2vr M 2 

2tt M 4 
6 



Im A 2 , 4 ,7,8(g 2 ,^) , 



(18) 



for the polarized structure functions. In the following we will consider the polarized part of T l 
only. 

For neutral current interactions the current operators obey 

j^ z] = j;> z . (19) 

Therefore, the crossing relation for the amplitude for q — > —q, P — > P yields 

TlM\-v)=TlM 2 ^) ■ (2°) 
The corresponding relations for the amplitudes Af c (q 2 , v) are 

= 4£ 8 ,9(9V), 
A^, 5 , 6 , 7 (q 2 , -u) = -A^ 5fiJ (q\v). (21) 

Furthermore, the amplitudes obey the following forward dispersion relations 

2 r°° , , v' 



^wtf.") = i/J/^l-^taV). (22) 

For the charged current interactions it is more suitable to study the linear combination of 
amplitudes 

TtM\ v) = T^(q 2 , v) ± T^ v \q\ v) . (23) 

Due to the transformation 

JT ] = JT > (24) 

the crossing relations read 

T ± (q 2 ,-u) = ±T ± (q 2 ,u) . (25) 
Correspondingly, one obtains for the combination of the amplitudes 

Af = Ai W ~ ± A t w+ (26) 

the relations 

^t 3 ,8,9(<? 2 , -v) = ±Af^8 :9 (q 2 ,v) , 

^4,5,6,7(9 2 , -V) = T^Wjl?^). (27) 

The respective dispersion relations for the amplitudes Af(q 2 , v) and A~(q 2 , v) are 

2 /■ 



v' 










V 








Z/2 




V 




V 12 




Z/2 




// 








I/ 2 



2 /"OO J/' 

^3,8, 9 (? 2 ^) = ~ / 2/ di/ ImA+ 3>8>9 (gV) , 

2 r°° 

2 f 00 i/ 
^,3,8,9(? 2 ^) = ~ J Q2/2 d"' u , 2 _ u2 lmA^ 8 9 (q 2 , v') , 

A-4,5,6,7(<7V) = ~J dl/^^lmA^^M. (28) 



In the case of charged current interactions we introduce the structure function combinations 

gt(x,Q 2 ) = gr-(x,Q 2 )±gf + (x,Q 2 ) . (29) 

The integral representations of the amplitudes Af c and Af can be finally expressed by the 
moments of the corresponding structure functions as 



and 



AM 2 



AM 4 



AM 2 



AM 4 



E 

n=0,2.. 



X n+1 JO 



dyy n g^ + (y,Q 2 ) , 



E ^I'dyy-g^Q 2 ), 



n=2A 



E 

n=l,3.. 

E 

n=l,3.. 



Jo 



JO 



dyy n g^ + (y,Q 2 ) , 



^i>, 9 (?V) 


AM 2 

V 


E 

71=1,3.. 


1 

j,n+l 


/ dyy n g 139 (y,Q 2 
Jo 




AM 4 

V 2 


E 

n=l,3.. 


1 

^,n+l 


1 dyy n g2 A (y,Q 2 ) 
Jo 




AM 2 

V 


E 

n=0,2.. 


1 

^n+l 


t dyy n g^(y,Q 2 ) 
Jo 




AM 4 

V 2 


E 

n=2,4.. 


1 

.£71+1 


1 dyy n g^(y,Q 2 ) , 
Jo 



(30) 



performing a Taylor expansion of Eqs. (|22| , [28|) and using Eqs. (|18| 
and the integration variable y = Q 2 / (2z/) . 



(31) 

where x = Q 2 /(2u) 



4 Operator product expansion 



The operator product expansion is one of the most general formalisms to analyze the properties 
of the structure functions in deep-inelastic scattering. We apply it to the T-product of two 
electroweak currents, 

f; u = T(J^(x)J?(0)) . (32) 

Near the light cone one obtains for neutral currents 

= q(xh fl (gv 1 +9A 1 l5)S(x)-f u (g V2 +g Aa -f 5 )P + q(0) 

+ q{^)lu{gv 2 + 9A 2 l5)S(x)^ ll (g Vl + g Al l5)P + q(x) , (33) 

and for the charged current combinations 

rp± rpW- I rpW + ( ^A} 

1 jlV ~ 1 fiU -'- 1 fill ' \°^) 



s 



where 



T. 



Q(x)jt,(g Vl + 9A 1 j5)S(x)-f u (g V2 + g Aa i5)P ± q(Q) 



± q(0hu(gv 2 + 9A 2 75)S(x)-f ll (g Vl + g Al -f 5 )P ± q(x) . (35) 
Here we used the projectors 



P + = 1, P~ = r 3 =\ . (36) 




The Fourier transforms of the expressions (|33]) and (|35|) have the following form 



d*xe iqx f^ c 



U{k)^ fl {g Vl + &Ai75) 7T + ^ o lu{9v2 + gA 2 ls)P + U{k) 



(2tt) 4 v ' ^ ,0/ (A; + g) 2 -m 

Tj^^h/^vi + ^Ts) ^ _ ™^ v(9Vz + g A2 j5)P + U (k), (37) 



and 



^ k - l U(k)-f ll (g Vl + g Al lh) n + m % 7v(^v 2 +^A 2 75)^ ± ^(^) 



( 2?r )4 V > '°'(jfe + g )2_ m 2 

T / (04 F (^)^(^ + SM^j^^^Yui^ + g M ^U{k). (38) 

Here m g denotes the mass of the outgoing quark in the scattering process and 

U(k) = J -0^ e-^(x) . (39) 

In the case of forward scattering the Fourier transform of the current product depends only on 
one variable. 

So far our discussion was quite general. In the following we will disregard the quark mass 
effects. As a consequence current conservation is also obtained for weak currents and only the 
structure functions Fi\^ =1 and #i| i=1 contribute. To isolate the spin dependent part in Eqs. ( p7|) 
and ( ]3"%D we use the identities 

= g^ag^p + g^gua — g^gap • (40) 

For the spin dependent part of T % one obtains 

TjL = - i (gv 1 gv 2 + g Al g A2 )e m upq a u1 + (g Vl g A2 + g Al gv 2 )S mU fi\q a ut + u af3 } , (41) 
T^u = ~ i (gv 1 gv 2 + g Ax g A2 )e m ufiq a vt + (g Vl g A2 + gA 1 gv 2 )S m vp[q a v+ + v al3 } , (42) 
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with 



ui = - 




u *p = _ 

J 


f i>k U(k) 


4 - -j 




= -j 





7/375 p+j^fc) 

P + C/(A;) - (g^-g) 



(£; + g) 2 



(k + q) 



(k + q) 



-P-U(k) T (q^-q) 
.p-U(k) - (?<--g) 



(43) 
(44) 
(45) 
(46) 



in the massless quark limit m q — > 0. Expanding the denominators (A; + g) 2 and (k — q) 2 in the 
above relations into powers of the ratio 2k.q/Q 2 the corresponding operators can be expressed 
in terms of the two operators B^lW'"'/*™} : 



7/< 8 



E ( 

n even, odd 



U 



a/3 



E 



/ 2 \ n+1 

2 \ n+l 



g w • • • g Mn e + ^-^> , 



4 



,a/3 



Q 5 



n even,odd 

2 \ n+1 



E 

n ocici 



where 



0±/?{>l-Mn} 



(2^ 



0-/9{a/ii-Ain} 



U(k) 1(3l5 k^...k fl ,P ± U(k) . 



(47) 

(48) 
(49) 

(50) 
(51) 



Inserting the expressions for the quark operators (|47H50l) into the spin-dependent part of the 
T-product, Tg pin , one obtains 



T+ • 

uv,spin 



-%Vi#v 2 + 9A 1 gA 2 )e fia upq a E 9 Ml ' ' ' ^" ( 7^ ) 

n even \W / 



n+l 



+ (9V 1 9A 2 +9V % 9Ai, 



-9r E 

n even / 

( 2 \ n+1 

n odd / 

( 2 \ n+1 

n ewen \ Q J 



0+Ml{M2— Mn} 



(52) 



and 



T, 



uv,spm 



<9 Vl 9v 2 + 9A x 9A 2 )e^q a Ef - f ^) Q-^"^ 

n odd / 

-9^ E <? Ml ■ ■ ■ - 

n odd \ 



+ (9v 1 9a 2 + 9v 2 9a 1 ) 



Q-"l{"2-«n} 
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/ 2 \n+i 

+ E ? w - • • ^ ™ (g M _ " {Ml "^ l} + g„e-^ 1 "^"> 

n ei>en \ / 

/ o \ n+1 

n odd \Q / 



(53) 



The operators B^^ 1 ""/M can be decomposed into a symmetric part, 0£, and a remainder, 0#, 
respectively 



} = Q±/3{M1-Mn} + Q±/3{/,l]-Mn} 



where 



6 



1 



e 



±/3{^i---^„} 
R 



n + 1 
1 

n + 1 



Q±/3{Ml"Mn} _j_ Q±Mi{/3 --Ain} _|_ _|_ Q±Mn{Ml-"/3} 



(54) 



(55) 



Q±(3{Hl---fl n } _ Q±Ml{/3 -'M„} _|_ Q±l3{m/12---/J.n} _ Q±M2{Ml/3-- Mn} _|_ ^gg^ 



In the massless quark limit both operators s and Or are traceless. To find the target mass 
dependence of the structure functions we have to construct the traceless nucleon matrix elements 
of these operators ||. 

The rank-n symmetric, traceless tensor constructed from n nucleon momenta P^ (P 2 = M 2 ) 
has the form 



n 



mi •■■/•»■ 



n ( p ) 



[n/2] 

E 

3=0 



2i n\ 



■ g...g P_^_PM 2j 

i n ~ 2 i 



(57) 



Here, [n/2] = n/2 for even n and [n/2] — (n— l)/2 for odd n. The sum contains j metric tensors 
g^Uf, with indices chosen in the set {ui, . . . , u n } in all possible combinations. The remaining 
n — 2j indices are carried by the momenta P^ r The sum contains n\/[2 J j\(n — 2j)\] terms. 

For the representation of the 0-operators traceless and symmetric tensors have to be con- 
structed from one nucleon spin vector 5* and n momentum vectors P. These tensors have the 
form @ 

M^-^(P ) S) = {S^P" 2 ...P^} 

= ^E^ nw '" [rt] "' Mn ( p ) - ^E^A nM1 ''' M ''' [Mj] '' M "( p ) > (58) 



where the indices in brackets [ui] indicate that the corresponding superscript is to be removed 
from the sequence. The coefficient of the second term in the above expression is chosen such 
that the trace vanishes identically. A more compact expression for these tensors was obtained 
in Ref. 101, 



M^-^P, S) = {S^P" 2 . . . P" n } 



f E 



n 



2> 



nl 



g---g [SP---P] S M 2 > .(59) 



n-2j 



Again the sum contains j times the metric tensors g^^- The remaining n — 2j indices are 
symmetrized in the product [SP . . . P]g. 
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5 The Nucleon Matrix Elements 



We calculate now the nucleon matrix element of the operators derived in the preceding section. 
The expectation value of the symmetric operators G5 is given by 



(ps\e^ { ^ n} \ps) 



(60) 



Correspondingly, for the operators 6# one obtains, cf. also 



n+l 



j^{"l«2'--«n} _ ^Ail{/3M2---Mn} _|_ 



j^-/3{«l "2 ■•■««} _ j^Mn{MlM2-'-/3} 



where 



(p, S) = ^Is^-^iP) + - Ar P^Mf-^fP, 5) . 
n+l n + 1 



(61) 



(62) 



In the massless nucleon limit the known results, cf. Ref. JT5| , are reproduced, 



{PS\Q^ { ^ n] \PS) 

(psio^^-^ips) 



n + l 



Ct/3 pHl pil2 pV-n _|_ pP pl^2 pV-n _|_ 



(63) 



n + l 



(gP p«l _ pP^j pM2 p^n 

+ (^S^P^ 2 — S^ 2 P^ p^ipV3 p^n 

_|_ _|_ fgP p^n _ g«n p!3^ p»i pP-i pUn-l 



(64) 



ct„ and d^ denote the twist-2 and twist— 3 operator matrix elements in the case of neutral current 
(+, NC) and charged current interactions (±) of moment index n. These are non-perturbative 
quantities and are independent of the nucleon mass. All information about the target mass 
corrections is contained in the tensors Mf' 11 '"^™ and M^ 1 '"^" etc. 

Let us consider first the nucleon matrix element of the first term of the operator T^,' spin , 
Eq. (R), 



++A 

fiv,spin 



Q 1 



(65) 



to derive the structure of the target mass correction for forward Compton scattering. We, 
furthermore, consider only the symmetric part of the operator 0. 



{PS I Pny spin 

\PS) 



uau/i 



[n/2] 



i ^ ^ x (n-2j+l) [ Q2 

[. n even j=Q ^ \ V 

[n/2] 1 

+ X/ ( n -2j+l) ( 7)2" 

n ewerc j=o \ ^ 



M 2 V' ^(n-j + l)! 



z/ j!(n-2j)! (n + l) 2 
M 2 \pP(S-q) (n-j + 1)! 



v 2 j\(n - 2j - 1)! (n + l) 2 



.(66) 
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Shifting the summation index from n to n = n + 2j both sums become independent and one 
obtains 



(PS\ T^' sp i n \PS) 



x 



-i{gv x gv-2 + 9A 1 9A 2 )e fia ui3q a 



2^ 2^ ^.fn+il 1 ra 



, ,p x^ n+1 ) \ Q 2 J v 2 j!(n-l)! (n + 2j + iy 



.(67) 



Carrying out the same steps also for the other parts of the expressions (^) and (|5^) and rear- 
ranging the tensor structures according to the structure functions one obtains for the twist-2 
part of the forward Compton amplitude the relations 



i n )j) a n+2j 



j=0 



q a P l3 (S ■ q) - q a S f3 (P ■ q) 



+ (gv 1 gA 2 + gv 2 9Ai. 



X 



j=0 



-g»v 



<mv \ s ■ g 

q 2 I v 



n odd X j=0 U ^ Z J 

S^P U + P^S V P^P V 



+ 



[(n + 2j + l)(n + 2j)+2j\ 



v 



(S-q) 



- E ^™E C (^ 

n odd j=0 



a n+2j 

n + 2j 



pa tju 9 1 00 

+ (S-q)- £ -^n(n+l)£C(n,j 



L n+2j 



V 



v " , x' L n + 2j 

n odd J=U ■> 



(6S 



and 



(PS\T ' ' in \PS) 



, q a S 13 1 00 

'^^2+ gA i gA 2 )\e i _ ia v(3— E ^Tl( n + ^E^^'i) fl n+2j 



V 



n odd 



X 



j=0 



q«pP{S-q)-q«SP{P-q) 1 « _ 



n odd 



i=o 



+ {gv^A 2 + gv 2 gA 1 )\ (-g^ + J V i )- 



E 



+ 



+ 



- £ C(n,j)^^[(n + 2j + l)(n + 2j) + 2j] 

-|-_L ' fl -} 9 

n even j=Q 

S»P V + P^S" P^P 1 



V 



-(S-q) 



4 



1 



- E 37™E C ( n ^ 



T 

n even j=0 



n + 2j 



(69) 



Correspondingly, the twist-3 contributions read 

(PS | T ui/,spin 

r=i \PS) = 



13 



( q a S^ 1 00 

i(9V 1 9V 2 + 9A i gA 2 )<£^aul3—— E Z^TT 4 E C ( n ' 3)3 d n+2j 

n even J=0 

g Q P^(^ • g) - g a ^(P • g) 



I/ 2 



E Z^T(^ + 1 )ECi(n-l,i)(n + 2i)d+ +2i 



x 



igv 1 gA 2 + gv 2 9A 1 )U-gau + ^)^- 1 E i 4 E c (^')j TT^rf^ 

L y U n odd- 1 j=0 K ,L ^^J) 



A — 2 — (5 - g) Jr; E ^ E 1( ,jKU 

n odd j=U 

x [(n - l)(n + 1) + 4j(n + 2j + 1)] 

p^ 



O 1 OO s 

(s-g)- E -E c i("- 1 >i)(^ + i + 1 )i rf n+2i , (70) 

V V ' n odd X j=0 ' 



and 



_ 3 f g Q D p 1 
(-PS I ^', S p7n I PS) = -i(g Vl g V2 + g Al g M ) < e mufi E "^I 4 E C ( n >3 )j d+ +2j 

L V n odd X j=0 

q a P (3 (S ■ q) - q a S?(P ■ q) 



1 °° 1 

E —^{n + l )Y. C ^ n -^3){n + 2j) d+ +2j \ 

n odd X j=0 } 



L n+2j 



00 



( S»P» + P»S» F"P" \2 1 

5 —(^-e) - E 37 E^- 1 ^)^- 



" n even j=Q 



x [(n - l)(n + 1) + 4j'(n + 2j + 1)] 

pfipv 



O 1 OO n 

(s-ff)- E -E c i("-U)(™ + j + i)j^ +2i , (7i) 

where the combinatoric factors C(n,j) and C\{n,j) are given by 

C( \ (n + j + l)! m , 

ln ' JJ W 2 / jln\(n + 2j + l)(n + 2j + 2y [ ' 6) 

In the presence of target mass effects all polarized structure functions receive twist-3 contribu- 
tions. Implicitly this effect was contained in Ref. |H| already in the case of the structure function 
9x- 
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6 Expressions for the Moments of Structure Functions 



To obtain the moments of the twist-2 and the twist-3 parts of the individual polarized structure 
functions we first consider the target mass corrections to the single amplitudes A^£l£'f. The 
twist-2 terms are : 



A?%(q 2 ,v) 



' i ■ ■ , j Q 



n even 
f2 



ml E M)!±M)!)„£ c( „,, )o+ , 



n+2j ? 



n even 
r2 Q„Q <l oo 



n odd 



t n+2j 



M 4 4<#<fl , ^vv, 

— E — ^ + i)Ec(M) n + 2j 



rt odd 



M 4 



1 „1 oo 



E ^EC(n,j)- : S[(n + 2j + l)(n + 2j) + 2j], (74) 



L n+2j 



n odd 



X 



3=0 



and 



Ai T=2 {q\v) = 

4r=2(?V) = 

A 3T=2 (q 2 ,v) = 
A^ T=2 (q 2 ,u) = 
A^ T=2 (q 2 ,u) = 

The twist-3 terms read : 

M 2 



M 2 ^ ((g q v ) 2 + (g q A ) 2 ) 



V 



E 



n odd X j=0 

( ( „q \2 i ( „q \2\ oo 



n+l 



n + 2j 



(n+l)J2C(n,j)a n l 2j , 



v 



n odd 



X 



3=0 



L n+2j 



M' 



^ E 



,2 

M 4 



3=0 



Lb n+2j 

n + 2j 



E %^EC(n,j)^fT[(n + 2 J + l)(n + 2 J ) + 2 J ]. (75) 



q oo 



i=0 



n + 2j 



A?%(q 2 ,v) 



A N 2 C Mq\») 



(q 2 , v) 



x 3 t=3 



M 
v 

M' 
v 



n euen *~ J=0 

2 tt„q\2 , / <?\2\ oo 



x i=o 



n even 

2 a „q „q oo 



E ^ E Ci(n - 1, j) [(n - l)(n + 1) + 4j(n + 2j + 1)] d&y , 

i odd X j=0 



n odd 



3=0 



n odd 



X 



j=0 



n + 2j 



(76) 
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and 

^2 t= 3 (<?V 

^=3(9 V 

A^ T=3 (q 2 ,u 



E U ^^J + r ; - 4£C(n,j)j^ 



n ocici 



X 



n+2j > 



M 2 _ ((^) 2 + (^) 2 ) 



i=o 

00 



27 jj £ 

ra oaa 



n+1 



■5](n + l)Ci(n-l ) j)(n + 2 J ')^ , 

j=0 



M 



M 



7v v 

T 

n enen J=0 

1 i«„9„8 



E E Ci(n - 1, j) [(n - l)(n + 1) + Aj(n + 2j + 1)] dfty 



,2 

M 4 



E + 1) E Ci(n - l,j)(n + j + l)jd% 2j 



j=0 



~ d+q 



E ^Ec(m) 



7 (ii ■ ■ rr*ll 

v n even ^ j=0 



n+2j 

n + 2j 



(77) 



From the representations Eqs. (^,^T|) the following relations between the moments of the 
twist-2 and twist-3 parts of the structure functions gf c,± (x, Q 2 ) and the operator matrix ele- 
ments a^ ,q and d^' q are obtained. The twist-2 contributions are 



f 1 dxx n g^ 2 (x,Q 2 ) = j: 
Jo „ 



(g q v) 2 + {g q A ? 
4 



in 



i)£c(n,j). 



(g q v) 2 + x 



E^,j) 



+q 



n+2j ) 



a 



fdx x n g 2 w %(x,Q 2 ) = -j: 

/ dx x n g^ 2 (x, Q 2 ) = 2g q v g« A (n + 1) £ C(n + 1, j) 

r l 00 

jf dxx n g^ 2 (x,Q 2 ) = Y: g q v g q A(n + l)(n + 2)J2C(n + l,j) 



C 2J , n=0,2... (78) 
ra = 2,4... (79) 



n+2j+l 



(ra + 2j + l) 



n = 0,2 ... (80) 



a 



n+2j+l 



/ cte x n g^(x, Q 2 ) = t;9v9a E C K ^7— ^ 
Jo „ 2 (n H 



i=o 

+9 



(n + 2j + 1) ' 

n = 2,4 ... 



-2j 



x[(n + 2j + l)(n + 2j) + 2j], n=l,3 



and 



/ (i«\ =2 (x,Q 2 ) = E 
Jo „ 



G#) 2 + {g q A ? 



in 



+ i)E c (^') 



n+2j ; 



(g q v) 2 + {g q A? x 



E c ( ra ^') ( 



/ x n £ 2r=2 (x, Q 2 ) = - £ 

-i 00 
y o ^ x^ 3 " r=2 (x, Q 2 ) = £ 2^(n + 1) £ C(n + 1, j) 



ra+2j ) 



n = 1, 3 



71 = 1, 3 ... 



(81) 

(82) 

(83) 
(84) 



a 



n+2j+l 



{n + 2j + 1) 



n= 1,3 ... (85) 



dx x^ 4 - =2 (x, Q 2 ) = £ <&fl(n + l)(n + 2) £ C(n + 1, j) " +2 ^ 



n = 1, 3 



(86) 



/' 
Jo 



1 00 a _<? 

x n g^ T=2 (x,Q 2 ) = £ £ C (n, j) ^ 



i=o 
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x[(n + 2j + l)(n + 2j)+2j], n = 0,2 
For the twist-3 contributions to the structure functions the moments are given by 



and 



[ l dx x n g?%(x,Q*) = Y: 
Jo „ 



(g q v) 2 + (g q A ) 



q\2 oo 



4 



■4£C(n,i)jdS 



+2j ) 



n = 0,2 .. 



.,9 oo 



/* dx x^l^x, Q 2 ) = £ ( ^ )2 ± ^ - l,j)(n + 2j)(n + l)dft« 

• /0 9 4 j=0 



n = 2,4 ... 



-1 oo 

/ dxx^l + 3 (x,Q 2 ) = E^E C iKj)^ + 2 J+ i[(n + 2j)(n + 2j + 2)+4j] 

J0 q j=0 



/ dx x n gt°=t{x, Q 2 ) = 4 E E C i ( n ' •?') ( n + 3 + 2 )jdn+2j+i 

Jo q j=0 



n = 0,2 ... 
n = 2,4 .. 



dx x"fe + (x,Q 2 ) = -2Ys9WaY,^3)3 



j=0 



d n +2j 
n + 2j 



[ dx x n g 1 r=3 (x, Q 2 ) = E 
Jo - 

[ dx x n g 2 ~ T=3 (x,Q 2 ) = E 
Jo - „ 



(g q v ) 2 + {g q A? v 



-Aj2C(nJ)jd~l 2j 



n = 1, 3 ... 



'/ ' 3=0 

q \2 i /„q\2 oo 



(9 q vY + (9aY 



Y i C 1 (n-l,j)(n + 2j)(n + l)d-l 2j , 



j=0 



n = 1, 3 



/ ^ x^ 3 - T=3 (x, Q 2 ) = E E Ci(n,i)d n+2i+1 [(n + 2j)(n + 2j + 2) + 4j] 



/ dx x n gl T=3 (x,Q 2 ) = 4^3^^Ci(n,j)(n + j + 2)jd n+2i+ i 

J0 9 j=0 



n = 1, 3 ... 

n = 1, 3 



dx x"# 5 r=3 (x,Q 2 ) = -2j29v9 q Aj2 C ( n ^)j 



j=0 



a n+2j 

n + 2j 



n = 0,2 



«0 



(89) 

(90) 
(91) 



n = l,3... (92) 



(93) 



(94) 

(95) 
(96) 

(97) 



Eqs. (|78|-|97|) can be used directly in structure function analyses based on integer moments. 
Unlike the massless case the operator expectation values and d^ do not decouple for a given 
spin n in this representation. As will be shown in the subsequent section, however, the infinite 
sums in Eqs. (|78|-|97|) can be related to integrals over one-dimensional partition functions via an 
analytic continuation from the integers n to the complex n-plane. Carlson's theorem PTJ assures 
the uniqueness of the analytic continuation and the inverse Mellin transform. 



7 The Inverse Mellin Transform 



In many practical applications, as the analysis of experimental data, the expressions for the 
moments of the twist-2 and twist-3 contributions to the deep inelastic scattering structure 
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functions at integer values of n are less suited than the corresponding x-space expressions. For 
this purpose we perform the inverse Mellin transform of the results derived in the last section 
and provide integral representations which can be directly applied to the measured structure 
functions, which are given as distributions in x and Q 2 , to unfold the target mass corrections. 

We apply a procedure analogous to that used by Georgi and Politzer in Ref. || in the case of 
unpolarized structure functions. We summarize the different steps for the case of the structure 
functions gf(x,Q 2 ). From Eqs. (|78| , |83D one obtains 



gf T=2 (x,Q 2 



E 



(g q v ) 2 + W 2 



M 



x— / dnx-( n+1 \n + l)J2 \Ta 



2\3 



[n + j + 1] 



l n+2j 



3=0 



Q 2 



]\n\ 



{n + 2j + If 



.(98) 



The operator expectation values are the moments of distribution functions F ±q (x), which are 

he massless li 

dyy n F^(y) . 



related to the polarized parton densities in the massless limit, Aq(x) ± Aq(x), cf. e.g. [15 

/' 

/o 



a 



(99) 



Hence, 



l n+2j 



(n + 2j + l] 



dyy n+ ^Gf(y) 



where 



Gf(y) 



E 



(g q v? + m 



q\2 -i 



h V\ hi V2 



(100) 



(101) 



Inserting these expressions into Eq. ( p8[ ) and interchanging the integrations and the summation, 
the following expression is obtained 



9lT=2( X ) 



/ dyGf(y)— dny-x-^in + lfY, 

Jo 2m J-ioo ~q 



My 

~Q 2 



2„,2\ J 



[n + j + l)\ 
j\{n + l)\ 



.(102) 



The infinite sum over j in this equation is performed by observing that 



(1 - r) n 



+i 



E 

Lk>0 



k 
I 



(103) 



This relation can be derived by differentiating the generating functional of the geometric series, 
k 



1 



r) . Here, 



denotes the Stirling numbers of the first kind, cf. |22 



The integration variable n is brought into exponential from by introducing differential oper- 
ators w.r.t. x 



9f T =2(x) 



d d 

doc doc 



dy 



Gf(y) 



x(l-M 2 y 2 /Q 2 )27ri J-ioc 



+ioo 



dny n x 



M l y 

~Q 2 



■ (104) 
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The n-integration leads to a 5-function 



Finally one obtains 



d d r 1 
x— x— 



dy- 



Gf(y) 



X 1 



M l y 



ln(y) — ln(x) — In 1 



2„,2 N 



M 2 y 

~Q 2 



(105) 



d d 
x—x— 

dx dx 



x 



Gt(Z) 



Here, £ denotes the Nachtmann-variable [|TTJ[] . 



(1 + 4M 2 x 2 /Q 2 y/ 2 £ 



2x 



l + (l + 4M 2 x 2 /Q 2 ) 1 /2 ' 
Similar expressions are obtained for the other structure functions 



03^=2 0*0 



05^=2 0*0 



ax z 



x 



(l + 4M 2 x 2 /Q 2 ) 1 / 2 £ 



2x' 



dx 2 



x 



Gf(0 



{l + 4M 2 x 2 /Q 2 y/ 2 £ 2 



(106) 

(107) 

(108) 
(109) 



2 d d 2 
ax dx z 



x 



—x 



dx 



(l+4M 2 x 2 /Q 2 y/ 2 £ 2 
Gf(C) 



(110) 



(l + 4M 2 x 2 /Q 2 ) 1 / 2 £ 



+ Q 2 X dx 2 



x 



Gt(0 



(l + 4M 2 x 2 /Q 2 ) 1 / 2 £ 



The functions G 2 (y) and G 3 (y) are related to the distribution function F ±q (y) by 

V J y J yi Vz J v* 



„ z jy yi 



(111) 

(112) 
(113) 



Performing the derivatives in Eqs. ( J106 , |108| - |111| ), the twist -2 contributions to the structure 
functions can be expressed in terms of the distribution functions F ±q (^) as follows : 



flT~T= 2 0) 



E 



+ 



4M 2 x 2 



q 2 e(i 

4M 2 x 2 (2 - 4MV/Q 2 



\£ (1 +4M 2 x 2 /Q 2 ) 3 / 2 
AM 2 x 2 /Q 2 ) 2 Ji Ci 



9ir=2{x) 



E 



Q 2 2(l + 4M 2 a; 2 /g 2 ) 5 / 2 y c & 4 6 ^ 2j _ 
G#) 2 + G&) 2 [ x ^(0 



(114) 



£ (1 + 4M 2 x 2 /Q 2 ) 3 / 2 



+ 



x(l-4M 2 i(/Q 2 ) f 1 ^! 



£(l + 4M 2 x 2 /Q 2 ) 2 Jt & 
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3 

+ 2(1 + 



4M 2 x 2 /Q 2 fd^ rid& F ± q( n ' 

4m%vq 2 )5/2 y 5 ^ 4 £ 2 



(115) 



9v9a' 



2x 2 



£(1 + 4M 2 x 2 /Q 2 ) 3 / 2 
6 4 6 



12M 2 x 3 /Q 2 f 1 <%! f 1 d& 



(l + AM 2 x 2 /Q 2 ) 2 
3 (4MV/Q 2 ) 2 
2(1 +4M 2 x 2 /Q 2 ) 5 / 2 

V oaf* 2 F^(0 



6 £2 £3 



1 ^3 F±g(6) - 



+ 



8M 2 x 2 



x(x + £) 



k £ (l + 4M 2 x 2 /Q 2 ) 2 Q 2 £(l + 4M 2 x 2 /Q 2 ) 3 
4M 2 x 3 (l-M 2 x(4x + Q/Q 2 ) figi f'd^ ± 
Q 2 (l + 4M 2 x 2 /Q 2 ) 3 Ji & 4 6 1 J 
3(4M 2 x 2 /Q 2 ) 2 (3-8M 2 x 2 /Q 2 ) /* f 1 ^2 f 1 d& F ± q{a ' 



(116) 

Si 



(1 + 4M 2 x 2 /Q 2 ) 7 / 2 



6 6 6 



1 



X 



+ 2 



+ 



f(l + 4M 2 x 2 /Q 2 

(i + O 



M 2 x 2 ( 
Q 2 l v £(l + 4M 2 x 2 /Q 2 ) 3 / 2 

£(l + 4M 2 x 2 /Q 2 ) 2 i ? 6 4 6 1?2J 

(1 + 4M%VQ 2 )V2^ ^ "Jf^ ^ F±9(6) )} • 

In the same way the moments of the corresponding twist-3 contributions, Eqs. 
inverted. One obtains the following x-space expressions : 



gf T=3 {x,Q 2 ) 
gf T=3 {x,Q 2 ) 
gf T=3 {x,Q 2 ) 

g^ T=3 (x,Q 2 ) 
gf T=3 {x,Q 2 ) 
Here we define 



4M 2 2 d 2 
~Q 2 ~ X dx 1 



x 



(l + 4M 2 x 2 /Q 2 ) 1 / 2 



x 



dx 2 



x 



(1 + 4M 2 x 2 /Q 2 ) 1 / 2 



d 3 



M 2 x 2 d 2 



x -— + 4— — — X-— -x 



dx 3 

M 2 3^ 

' A ~Q 2 ~ X dx^ 

M 2 2 d 2 

-9 T 

Q 2 dx 2 



X 



^(0 



Q 2 dx 2 J L(l + 4M 2 x 2 /g 2 )V 2 £ 



X 



(l + 4M 2 x 2 /Q 2 )V 2 £ 



x 



^(0 



(l + 4M 2 x 2 /g 2 )V 2 f 



G#) 2 + (0. 



<n2 h 



yi •% 1/2 



o ^ 3/ "'i/i 2/2 ^1/2 2/3 



(117) 



(118) 
are 

(119) 
(120) 
(121) 
(122) 
(123) 



(124) 
(125) 
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The matrix elements of the twist-3 operators are the moments of the distribution function 
D ±q (x) in the massless limit, 

d±i = ( l dxx n D ±q (x) , (126) 



which has, however, no partonic interpretation. Performing the derivatives in Eqs. ( |119| - |123| ) the 
following representation for the twist-3 parts of the polarized structure functions are obtained : 



4 Q 2 { (1 + 4M 2 x 2 /Q 2 ) 3 / 2 

3 fl ^(6) 



(l + 4M 2 x 2 /Q 2 ) 2 Jt 6 
(2-AM 2 x 2 /Q 2 ) f^dCi f 1 d^2 



{i + AM 2 x 2 /Q 2 fi 2 J^ 6 4 6 £>±9(6,) J ' (127) 



4 [ (1 + AM 2 x 2 /Q 2 f/ 2 

2J2 ln2 



1 - 8M 2 x 2 /Q 2 r'dti 
(l + 4M 2 x 2 /Q 2 ) 2 Jt 6 

- {i + AM 2 x 2 /Q 2 fi 2 h tL~^ (6) r ( 8) 

a± (*Q 2 ) - Wq4 2x(l -MV/Q 2 ) fdh^, 



Q 2 (l + 4M 2 x 2 /Q 2 ) 2 7 ? £i 4 £ 2 1 ^ 
,M 2 x 2 (1-4MV/Q 2 ) /-i ridfj /- 1 d£ 3 



+ /< «■ i f / 6 f b±,&) ) • (i29 » 



9 

„2 



Q 2 1 (l + 4M 2 x 2 /Q 2 ) 2 



x 2 (5 - 2m 2 </q 2 ) r 1 ^ ± 

(l + 4M 2 a; 2 /Q 2 ) 5 / 2 ./ 5 6 

6a; 2 (2x 2 - 3£ 3 ) r 1 f 1 d& n ± q(n 

e{l + ^M 2 x 2 /Q 2 fJi 6 4 6 1 J 

M 2 x 2 (l-MV/Q 2 ) /- 1 d£ 2 /-i d£ 3 



gir -M - -2S^^{ (1 + 4 ^ vg2)3 4 'f^(6) 

2(1-MM/Q 2 ) <"j6 ± „, r , 

6MW ^jW^fcjl USD 



(1 + 4M 2 iV0 2 ) 5/2 ' / 5 4 6 4 & 
The multiple integrals in the above expressions may be simplified further to single integrals by 
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using the idendities 



e 6 6 



0(6) 



■k -h: 6 ^6 43 A 6 



d-e-eiog^l 



0(6) • 



(132) 
(133) 



We would like to comment on the behavior of Eqs. (|1 14HTT8|J127HT31~D in the range of large 
values of x. At first sight these relations for the structure functions are not correct. The structure 
functions should vanish at x = 1, which is not the case for the r.h.s. of Eqs. (|114| - pT8|J127| - [r3TD . 
The inconsistency of these equations can be shown also by simply resumming the (M 2 / Q 2 )-terms 
in Eqs. (|78^p7|), see also p3| . As an example, for the moments of gf(x, Q 2 ) one obtains 



f 1 dxx n gt(x,Q 2 ) = [ 
Jo Jo 



i/(i-MVQ 2 ) d d 

dxx -r^— 
o ax dx 



x 



Gt(0 



[l + 4M 2 x 2 /Q 2 y/ 2 £ 



(134) 



The integrand in the r.h.s. of Eq. (|134|) is just the expression of gf(x,Q 2 ) after the inverse 
Mellin transform, Eq. ( |1 14|) . Therefore, the x-space expression for gi(x,Q 2 ), Eq. ( |114j) , cannot 
be completely correct because of the different ranges of integration. 

This problem can be resolved if we assume that F(£), and consequently Gi(£), vanishes for 
£ > Cthi where 6ft — 6 X = !)• This assumption is equivalent to the introduction of a kinematic 



threshold factor jnj in the moments of distribution functions, as was proposed in Ref. p4" 
Moreover, only in this case both the methods of Refs. 



1C1 are consistent. The matrix elements 



of operators calculated by Nachtmann moments of the structure functions ( |114j - |ll8D coincide 
with the definition of the matrix elements through the moments of the function F(x), Eq. (j99|) , 
as was observed for the unpolarized structure functions in Ref. [25]. In the range of large values 
of x the contributions due to dynamical higher twist operators are increasingly important. As 
well-known the number of higher twist operators grows with the spin index n the more the 
larger the twist r Q and contribute significantly because of their number as x — > 1. Therefore the 
moments of the distribution function F(x) cannot be expressed by the matrix elements of the 
lowest twist operators onlyQ. Conversely, the matrix element of the operators in the expressions 
of the moments, Eqs. (f78|- |87l) , cannot be approximated by Eq. (|99| ) alone. As a consequence it 
turns out, that in the approximation considered the results for the structure functions in x-space 
are reliable only at £ <^ ^ th . 

After having performed the resummation of the M 2 /Q 2 effects in Eqs. (114-TT£ , 127| -131) 
one may try a perturbative treatment if M 2 /Q 2 <C 1 by expanding these expressions into Taylor 
series. Let us consider the principle relations for the structure function gi(x, Q 2 ) as an example. 
Its moments are given by 



M[ 9l ] (N,Q 2 



o 



dxx n gi(x, 



(135) 



3 The number of these operators N(n, r) is a complicated function in most field theories, cf. e.g. || 0. 

4 Quite a variety of empirical ansatze were tried in the literature to model higher twist contributions to structure 
functions, see [p6[ p7[ and references therein. They are strongly correlated to the parameters of the lowest twist 
distributions |27| . These terms are not derived in QCD. Due to the structure of higher twist operators and 
coefficient functions, cf. Refs. [5-7], it is expected that the structure of the higher twist contributions is even 
more complicated than assumed in the ansatze quoted above. 
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The first terms of the corresponding Taylor series are 



9i 



+ 



+ 



#10+ q 2 

m 3 \ 



M 2 ( (n+l)(n + 2)(n + 5) n 



,n+ 1 



(n + 3) 2 



„n+2 I A „n+2 

#10 "T ^" 



n + 3 



#20 



(n + l)(n + 2)(n + 3)(n + 9) +4 (n + l)(n + 2)(n + 3) ^ 

fio / ; 7T7~ ; c \ #20 



2(n + 5) 2 



(n + 4)(n + 5) 



(n + l)(n + 2)(n + 3)(n + 4)(n + 13) 6 
6(n + 7) 2 510 

| ^ (n+ i) (n + 2 )(n + 3)(n + 4) n+6 - 
+ 2(n + 6)(n + 7) feo , 



0(f),d36) 



Here g>" and g^o are the moments of the corresponding structure functions in the limit M — > 0, 
with 



#io 

#2^ 



n 

n + 1 



(d n - On) . 



Let us consider furthermore the twist-2 part of g\{x) only,Q 



9i 



„ M 2 f(n + l) 2 (n + 2)^ +2 



#io + 



Q 2 



'M 2 \ 2 



(n + 3) 2 

\2/ 



#10 



+ 



,q 2 

'M 



2\ 3 



(n+l) 2 (n + 2)(n + 3) ^ n+4 

2(n + 5) 2 ^ 10 
( n + l)2 (n + 2)(n + 3 )(n + 4) 



6(n + 7) 2 



The inverse Mellin transforms of the prefactors of giQ~ 2 and (7i +4 in Eq. ( |138| ) are 



M- 1 



M 



-i 



(n + l) 2 (n + 2) 
(n + 3) 2 
(n + l) 2 (n + 2)(n + 3) 
(n + 5) 2 



(137) 



(138) 



(x) = - \x 3 — + 5x 2 J 5(1 - x) + Ax 2 In x + 8x 2 , (139) 

/ d 2 d \ 

f x ) = [x 5 — + 12x 5 — + 64x 4 U(l-x) 
\ dx 2 dx y 

-96x 4 lnx - 128x 4 , etc. (140) 



One observes that the higher the exponent k of (M 2 /Q 2 ) k the higher the order of the derivatives 
of the 5-function in the inverse Mellin transforms given above. Let us assume that the large x 
behavior of gio(x) is typically being described by 

g w (x) oc (1 - x) a , (141) 

with a > 0. The target mass effects, if represented in terms of a Taylor expansion, yield 



gi(x,Q 2 



1\/T- 

1 - x) a + — {ax(l - x) a ~ l - 5(1 - x) a + h(x)} 



Q 



(142) 



+ 



m 4 n 



Q 4 12 



a(a - l)x 2 (l - x) a " 2 - 7ax(l - xf + 31(1 - x) a + J 2 (x) \ + ... 



3 The subsequent arguments are similar in the case of the twist-3 contributions. 
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The integrals 



h(x) 



h(x) 



4/;^ [2+ ,„ (Il) ]( 1 -± 

-16 [4 + 31ii(n)] C 1 - — V 

Jx X\ \ X\) 



(143) 



are regular. However, starting at O(K) with K > a contributions oc 1/(1 — x)\ K a ' emerge, 
which diverge as x — > 1. The resummed expressions, on the other hand, are convergent. 

8 Relations between the Structure Functions : Twist 2 



In previous analyses three independent relations between the twist-2 parts of the polarized struc- 
ture functions were derived in the lowest order of the coupling constant, the Dicus-relation 



the Wandzura-Wilczek relation | 17j , and a third relation given in Ref. [TE\. We investigate in 
the following the impact of the target mass corrections on these relations. 

The expressions for the twist-2 contributions to the structure functions gi(x, Q 2 ) and #2(^5 Q 2 )-, 
Eqs. (11061 ) |108|) , can be rewritten in the following form 



gi(x,Q 2 ) 
92(x,Q 2 ) 



d d 2 
x—T(x,Q 2 ) + x 2 — -zFix, Q 2 ) , 
dx dx z 

-2x^F{x, Q 2 ) - x 2 -^F(x, Q 2 ) 



Here T{x, Q 2 ) is given by 



X 



'1 + 



AM 2 x 2 C 

~Q 2 ~ 



From Eqs. (|144 , |145| ) one obtains 

92(x, Q 2 ) 

and 

d 



-gi(x,Q 2 ) - x^-T{x,Q 2 ) 
dx 



x^-F(x, Q 2 ) 
dx 



Jx y 



from which the Wandzura-Wilczek relation [17[ 

92(x,Q 2 ) = - 9l (x,Q 2 ) + f l ^ gi {y,Q 2 

Jx y 



(144) 
(145) 

(146) 



(147) 



(148) 



(149) 



results. Therefore the Wandzura-Wilczek relation holds in the presence of target mass correc- 
tions. This result was obtained in Ref. before. The Q 2 -behavior in Eq. (|149|) refers to the 

I) one obtains 



(M /Q )-corrections in all orders. From Eqs. (|78l , [791, p3| 



n [ dxx n g{(x,Q 2 
Jo 



-in + 1) / dxx n g l 2 (x,Q 2 
Jo 



(150) 
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where n — 2, 4... for i = NC, + and n — 1, 3... for i — — . Note that the moments in Eq. ( |1 50|) 
are Cornwall-Norton moments, [E9]. Moreover, it can be shown that the Wandzura-Wilczek 
relation is also valid for the quarkonic operators even in the massive quark case, see appendix. 
The expression for g2, T =2{ x ,Q 2 ), E<1- (|15U|), is formally consistent with the Burkhardt- 



Cottingham sum rule [29] in the presence of target mass corrections 



dxg\ 2 (x,Q 2 ) = 



;i5i) 



although the Oth moment of the structure functions g 2 (x,Q 2 ) is not described by the local 
operator product expansion. The same holds for the Oth moment of the twist-3 contribution to 
g\(x,Q 2 ), cf. Eqs. OH). 

The relation between the twist— 2 contribution to the structure functions g${x) and g±{x) 
obtained in Ref. \TE\ in the massless limit is also valid at any order of M 2 /Q 2 . To show this we 
integrate Eq. (|110|) which results into 



gl(x J Q 2 )=2x f^g^y^Q 2 ) 

Jx y 



(152) 



This relation can also be derived directly from the relation between the moments of the twist-2 
parts of the structure functions g^{x) and g^x) 



n + 2 



dxx n g\{x, Q 2 ) 



(153) 



where n = 2, 4... for i = NC, + and n — 1, 3... for the i — — . 

In the limit M -> Eqs. (0, [S2|, |S§, 0) result into the Dicus-relation 



dxx n g l 4 (x,Q 2 ) = 2 I dxx n+1 g\{x,Q 



cf. []T6|, where n = 2, 4... for i = NC, + and n = 1, 3... for i 

gl(x,Q 2 ) = 2xgl(x,Q 2 ) 



-, which reads in x-space 



(154) 



(155) 



This relation is not preserved in the presence of target mass corrections as in the corresponding 
case for unpolarized deep inelastic scattering also the Callan-Gross relation [^DJ. This is expected 
since the tensor-structure in W^ v in the case of g§ and g^ is the same as that of F\ and F 2 for 
the unpolarized part except the factor of S ■ P. The breaking term Ad(x, M 2 /Q 2 ), 



M 2 



gl(x, Q 2 ) = 2xgl(x, Q 2 ) + A D lx, 



(156) 



is of 0(M 2 /Q 2 ) and reads 



M 2 

~Q 2 



2 d d 2 

M 2 , d 2 

—2 r 3 

Q 2 dx 2 

9 d 
+2x 2 — 
dx 



x 



[l + 4M 2 x 2 /Q 2 y/ 2 e 



x 



{l + AM 2 x 2 /Q 2 y/ 2 £ 
[l + 4M 2 x 2 /Q 2 y/ 2 £ 



(157) 
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9 Relations between the Structure Functions : Twist 3 



9.1 General Relations 



In the presence of target mass corrections all structure functions gi\ i=1 contain twist-3 contri- 
butions. On the contrary, in the massless limit this is the case for the structure functions and 
g 3 only |0|. The structure functions are defined by the hadronic tensor W^, Eq. (H). From its 
structure alone it cannot be concluded which structure functions vanish in the limit M — > 0. The 
mass dependence of the scattering cross sections for longitudinal nucleon polarization, Eq. (|T0|), 
however, reveals that both the structure functions #2 and g 3 do only contribute at 0(M 2 /Q 2 ). 
Therefore, a complete account for twist-3 contributions requires to consider also the nucleon 
mass corrections for the other polarized structure functions. 

Let us consider now the relations between the twist-3 parts of the structure functions. From 
Eqs. ( |119| - |123"D one derives the following new relations between twist-3 parts of different spin- 
dependent structure functions : 



4MV 

~Q 2 



AM 2 x 2 



Q 2 



g{, T =3^,Q 2 ) 

-9\, r=3(x,Q 2 ) 



A r= 3 (^Q 2 )-2 f -92. r=s(y,Q 2 ) 

Jx y 



9i r. 3 (*,Q 2 ) 1 + 



Jx y 



Q 2 



2xgi T=3 {x,Q< 



1 d y i ( n 2, 

x y 



(158) 
(159) 
(160) 



Here, Eq. ( |158| ) at the one side and Eqs. ( |159| , |16CI| ) on the other side correspond to different flavor 
combinations among the twist-3 contributions. This is similar to the case of the twist-2 terms, 
where the former case corresponds to the combination Aq + Aq, and the latter to Aq — Aq. 
The corresponding relations for the Mellin moments of the twist-3 part of the structure 



functions can be derived from Eqs. (p8H97|) and read 



dx x n gl T=3 (x,Q 2 



1 AM 2 x 2 

dx xn —Fy2—3l. T=3 (^<2 2 
Q z 



dx x + g§ T _ 3 (x, Q 



n + 1 /-i n 4M 2 x 2 i , _ 2 , 

dX X Q2 T=s( X ,Q ) 



n + 3 Jo 
n + A n 



(161) 



n + 1 Jo 



dx x n g\ T=3 (x,Q 2 



/•1 AM 2 r 2 
+ j Q dx x n —^9\ T =3^Q 2 ) , (162) 



doc x cj^ T= ^(yOc ; 



(163) 



2(n + 1) Jo 

The consistency of these relations is easily checked by partial integration. Eqs. (|158| - |160"D show 
that the twist-3 contributions to gi,g4 and g$ vanish in the limit M — > 0. On the other hand, 
if one keeps terms of 0((M 2 /Q 2 ) ■ #2(3)) and twist-3 in the scattering cross sections, one has to 
account also for the twist-3 terms in gi,g$ and g 5 . Because of the arguments given at the end 
of section 7 this requires the complete account for target mass corrections. 



9.2 Other Relations 

If the nucleon mass effects are disregarded in the operator product expansion only the structure 



functions #2 and g 3 receive twist-3 contributions [15]. A relation between the matrix elements 
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d^~, d q n contributing to g 2 and respectively, can be constructed by assuming that 



d q ~ = d q+ 



val 



One obtains 

dxx r 



o 



4# 5 (x>Q 2 ) - ^^93 {x,Q 2 ) 



X 



(n-1) 



dxx n 



ngJ{x,Q 2 ) + (n + l)g^x,Q 2 



vn—vp 



■yp—-yn 



ju- _ Jd- 



(164) 



(165) 



n 
12 



< + -< + ] val , n = 2,4,... (166) 



Eq. ( |164| ) allows to combine these relations to 



gZ=7(x,Q 2 ) = l2 



x92,t=3(x,Q 2 ) - I dyg 2 , T =3(y,Q' i 

J X 



jp—-yn 



(167) 



The target mass corrections break this relation. The correction terms to Eqs. (|165yi66|) are 



/ dxx n 
Jo 



4g 5 (x, Q 2 



n + 1 



x 



9 3 ( x , 



oo / M 2\3 
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(n + j)\ d 9 n+2j 



^\Q 2 J j\(n-l)\n + 2j 



doCtJc 



ngJ(x,Q 2 ) + (n + l)g](x,Q 2 



[n + j)\(n + 2j) C 2j 
j\{n-l)\ n + 2j 



:i68) 



. (169) 



In general it is not possible to absorb these contributions into the respective structure functions. 
A sum rule for the first moment of the valence part of gi(x, Q 2 ) and g 2 (x, Q 2 ) 



gX{x) + 2gX{x) 



0. 



(170) 



was obtained in Ref. |31 |. In the zero-mass limit it was shown in Ref. |15j that this sum rule 
is formally consistent with the operator product expansion. The valence parts gY(x) and g^ix) 
cannot be isolated for electromagnetic interactions from the complete structure functions. One 
may, however, refer to the charged current case, Eqs. (|83U84] ), and obtains 



f dxx{g 1 (x,xQ 2 ) + 2g 2 (x,Q 2 )) = V 
Jo " „ 



(9 Q v) 2 + (9 Q a) 2 ^ ( M 2 ^ 

4 hW 



(j + iK7 +1 . (i7i) 



The l.h.s. of Eq. ( |1 71| ) contains only valence quark contributions. For individual quark flavors 
we may write separately 



dxx 



gl«(x,xQ 2 ) + 2g?(x,Q' 



M 



2\3 



j=0 



(172) 



The matrix elements of the twist-3 operators for odd indices, d 2 j + \, are proportional to the mass 
of the quark. Let us consider the matrix elements of the following bi-local operators 



Qpm-nn = S n qvys<T P,t HD>* . . . iD^q , 



(173) 
(174) 
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where the symbol S n denotes the symmetrization of the indices fii,.. 
elements of these operators are related by 

(PSIQ^-^IPS) = m q (PS\0^-^\PS) 



,fji n . The nucleon matrix 



(175) 



In the massive quark case the 0- and O-operators are not traceless and their nucleon matrix 
elements are defined as follows 

(PS\Q^-^ \PS) = S n d 2 q J+1 {S?P^ . . . - S^PP . . . PW) + g l3 terms , (176) 
(PS\O p,xl -^\PS) = S n ^^-(S^P fil . ..P^ - S^P P ...P^) + g tj terms . (177) 

The matrix element in the r.h.s. of Eq. ( |177p is related to the moments of the twist-2 structure 
function h\(x), cf. Ref. ]32| , by 



dxx n hi(x) 



n-l 



hi{x) 



\n-l 



h\(x) 



(178) 



Due to the same symmetry properties of operators and O the tensors in the r.h.s. of the 
Eqs. ( |1 76| ) and (177) are the same. Therefore, using the expressions for nuclear matrix elements of 
the above operators and the relation between them, Eq. (172), the first moment of gx(x)+2g 2 (x), 
Eq. ( |171 ), can be expressed through the moments of transversity distribution 



dx x [gi{x] + 2g 2 (x)] 
After resummation one obtains 



£q Tflq 



2 M 



E 

3=0 



'M_ 

M 1 



2\ 3 



(j + 1) / dx x^lhix) - hxix)} . (179) 



/ dxx \qi(x) + 2q 2 (x)] = — — - / dx- 
Jo [yu ; y v n 2 M Jo 



h(x) — h(x) 
T M 2 x 2 \ 

K Q 2 



180) 



For Q 2 > M 2 , which was implicitly assumed in performing the light cone expansion above, the 
integral in the r.h.s. of Eq. ( |180| ) is finite. Therefore the r.h.s. of Eq. ( |180| ) vanishes in the limit 
m q — > 0. This is equivalent to d 2 j + i — > and Eq. (|170|) is found to be formally consistent with 
the the result of the operator product expansion. 



10 Conclusions 

We have calculated the target mass corrections for all polarized structure functions for both neu- 
tral and charged current deep inelastic scattering in the case of conserved currents. The results 
were obtained by using the local light cone expansion of the Compton amplitude for forward 
scattering. In the evaluation we used the approach of Ref. 0. The target mass corrections 
imply besides the twist-2 terms twist-3 contributions for all polarized structure functions. Only 
after the inclusion of the target mass corrections the description of the polarized scattering cross 
sections can be regarded as being completed at the twist-3 level, since the contributions to gi, g 4 
and g 5 are of the same order as those to g 2 and g 3 discussed previously. The corrections were 
both represented in terms of the integer moments which result from the light cone expansion 
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and their analytic continuation and Mellin inversion to x-space. The latter representations are 
directly applicable in experimental analyses. 

We investigated the effect of the target mass corrections on the sum rules connecting the 
polarized structure functions in lowest order in the coupling constant. For the twist-2 contribu- 
tions both the Wandzura-Wilczek relation |17 and the relation derived in Ref. |15| are preserved, 
whereas the Dicus relation \ Wj , similarly to the Callan-Gross relation |3(J in the unpolarized 



case, receives a correction. It was also shown that the Wandzura-Wilczek relation is preserved 
in the presence of quark-mass corrections. A previously derived integral relation between the 
twist-3 valence contributions to g 2 and #3 [H| receives target mass corrections. A relation for 
the first moment of a combination between gi and g 2 |U| is preserved. 

Three new integral relations were derived for the twist-3 contributions of the polarized struc- 
ture functions. They hold without further assumption on the flavor combinations of the related 
structure functions. In the case of the relation between the twist-3 contributions to g\ and g 2 
experimental test can be performed in the foreseeable future by precise measurements of the lon- 
gitudinal and transversely polarized deep inelastic eA-scattering cross sections at lower values 
ofQ 2 . 
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11 Appendix : 

The quark mass contribution to the Wandzura-Wilczek 
relation 

We consider the parts of forward Compton scattering amplitude T l spin , Eqs. (0) and (|53|), 
which correspond to the structure functions gi(x) and #2 (2) in the massive quarks case|] 

ft ~ B^q a E g^...g^(— ^_^) n+1 ± ^-^> . (181) 

n even. odd 1 A • q 

Here mx is the mass of the struck quark and m q the final quark mass. The initial-state quark 
is not necessarily on shell, but we assume that its wave function is the usual free-particle Dirac 
spinor for a particle with mass mx- 

In the massive quark case the 0-operators are not traceless anymore and have no definite 
twist. To construct the nucleon matrix elements for these operators, we have to decompose them 
into traceless ones. The most general form for a such decomposition is 



oo 



e ±/3{ m ... M „} = Y / A(n,j)0 ±t3{fll - fln ^ } g...g 



(m 



3=0 

3 



j ~° 3 

oo 

+ J2 C ( n >j) 0± ^ V '^- 2j} 9^9^A m x) j ■ ( 182 ) 

Jf ~° 3 

Because of the antisymmetric tensor e^avp the third term oc C(n,j) does not contribute to Tt- 
Symmetrizing the other two operators in Eq. (|182fl we separate the twist-2 part in the forward 



Compton scattering amplitude. The symmetric part of the remaining two operators are identical. 
Therefore one obtains 

[n/2] / 2 \ n+1 

ft ~ E E L m 2 —2 ) [A(n,j) + B(n,j)} 

n even,odd j=0 \ V m X + m q ) 

xf ■■• q »nQ±P{ f *i-Hn-v} g g{m\y 

3 

X [A(n,j) + B{U,J)} q" 1 ■ . . g Mn-2 J@ ±/3{«l -«n- 2j } 

4:m 2 x Q 



1 i elet,odd\Q 2 - m x + rn?J Xl 



'+1 r a 2 r>2 



(Q 2 — m\ + m 2 ) 2 



X 



q n . . .^e%i-w} . (183) 



In the last step we changed the summation index from n to Z = n — 2j and introduced the 
function Xi(z), 

X l {z) = J2(- Z y[A(l + 2j,j) + B(l + 2j,j)\. (184) 
3=0 



The resummation of quark mass effects was considered for the unpolarized structure functions in Refs. fa, B3| 
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As can be seen from Eq. ( |183|) we obtain the same tensor structure as in Eq. (|65|) which leads 
to the Wandzura-Wilczek relation. Because the dependence on the quark masses is the same in 
the structure functions g±(x) and gzi^x) the Wandzura-Wilczek relation is not violated by these 
corrections. 



31 



References 



[1] K.G. Wilson, Phys. Rev. 179 (1969) 1699; 

R.A. Brandt and G. Preparata, Fortschr. Phys. 18 (1970) 249; 

W. Zimmermann, in: Elementary Particle Physics and Quantum Field Theory, Brandeis Sum- 
mer Inst., Vol. 1, (MIT Press, Cambridge, 1970), p. 397; 
Y. Frishman, Ann. Phys. 66 (1971) 373; Phys. Rep. C13 (1974) 1. 

[2] A. Bodek et al, Phys. Lett. 30 (1973) 1084; Phys. Rev. D21 (1979) 1471; 
J.S. Poucher et al., Phys Rev. Lett. 32 (1974) 118; 
E.M. Riordan et al., Phys. Rev. Lett. 33 (1974) 561; 
S. Stein et al., Phys. Rev. Dll (1975) 1884; 
W.B. Atwood et al., Phys. Lett. B64 (1976) 479; 
M.D. Mestayer et al., Phys. Rev. D27 (1983) 285; 
R.G. Arnold et al., Phys. Rev. Lett. 52 (1984) 727; 
S. Dasu et al., Phys. Rev. Lett. 61 (1988) 1061; 
R.E. Taylor, Rev. Mod Phys. 63 (1991) 573; 
H.W. Kendall, Rev. Mod. Phys. 63 (1991) 597; 
J.I. Friedmann, Rev. Mod. Phys. 63 (1991) 615; 

L.W. Whitlow et al., Phys. Lett. B282 (1992) 475; SLAC-report-357 (1990) and references 
therein; 

W. Bartel et al., Phys. Lett. 28B (1968) 148; 
W. Albrecht et al., Nucl. Phys. B13 (1969) 1; 
J. Moritz et al., Nucl. Phys. B41 (1972) 336; 

E. Lohrmann, in : Proc. Lund International Conf. on Elementary Particles, ed. G. von 
Dardel, Lund, Sweden, June 1969, (Berlingska Bocktryckeriet, 1969), p. 13. 

[3] D. Adams et al., SMC collaboration, Phys. Rev. D56 (1997) 5330; 

B. Adeva et al., SMC collaboration, Phys. Rev. D58 (1998) 112001; 

K. Ackerstaff et al., HERMES collaboration, Phys. Lett. B404 (1997) 383; 
A. Airapetian et al, HERMES collaboration, piep-ex/9807015| . 

[4] K. Abe et al., E143 collaboration, Phys. Rev. Lett. 76 (1996) 587; Phys. Rev. D58 (1998) 
112003; 

K. Abe et al, E154 collaboration, Phys. Lett. B404 (1997) 377. 

[5] S. Gottlieb, Nucl. Phys. B139 (1978) 125; PhD thesis A New Twist in Deep Inelastic 
Scattering, Princeton University, 1978. 

[6] S. Wada, Progr. Theor. Phys. 62 (1979) 475; Nucl. Phys. B202 (1983) 201; 
M. Okawa, Nucl. Phys. B172 (1980) 481; B187 (1981) 71; 
H.D. Politzer, Nucl. Phys. B172 (1980) 349; 

E.V. Shuryak and A.I. Vainstein, Phys. Lett. B105 (1981) 65; Nucl. Phys. B199 (1982) 
451; B201 (1982) 141; 

R.L. Jaffe and M. Soldate, Phys. Lett. B105 (1981) 467; Phys. Rev. D26 (1982) 49; 
S.P Lutterell, S. Wada, and B.R. Webber, Nucl. Phys. B188 (1981) 219; 
S.P Lutterell and S. Wada, Nucl. Phys. B197 (1982) 290; 

R.K. Ellis, W. Furmanski, and R. Petronzio, Nucl. Phys. B207 (1982) 1; B212 (1983) 29; 

C. S. Lam and M.A. Walton, Can. J. Phys. 63 (1985) 1042; 



32 



A.P. Bukhvostov and G.V. Frolov, Sov. J. Nucl. Phys. 45 (1987) 704; 
J.W. Qui, Phys. Rev. D42 (1990) 30. 

[7] J. Bliimlein and W.L. van Neerven, to appear. 

[8] H. Georgi and H.D. Politzer, Phys. Rev. D14 (1976) 1829. 

[9] A. Piccione and G. Ridolfi, Nucl. Phys. B513 (1998) 301. 

[10] O. Nachtmann, Nucl. Phys. B63 (1973) 237. 

[11] V. Baluni and E. Eichten, Phys. Rev. Lett. 37 (1976) 1181; Phys. Rev. D14 (1976) 3045; 
A. De Rujula, H. Georgi, and H.D. Politzer, Ann. Phys. (NY) 103 (1977) 315; 

D. J. Gross, S.B. Treiman, and F.A. Wilczek, Phys. Rev. D15 (1977) 2486; 
A. De Rujula, H. Georgi, and H.D. Politzer, Phys. Rev. D15 (1977) 2495. 

[12] S. Wandzura, Nucl. Phys. B122 (1977) 412. 

[13] S. Matsuda and T. Uematsu, Nucl. Phys. B168 (1980) 181. 

[14] H. Kawamura and T. Uematsu, Phys. Lett. B343 (1995) 346. 

[15] J. Bliimlein and N. Kochelev, Phys. Lett. B381 (1996) 296; Nucl. Phys. B498 (1997) 285. 

[16] D.A. Dicus, Phys. Rev. D5 (1972) 1367. 

[17] S. Wandzura and F. Wilczek, Phys. Lett. B72 (1977) 195. 

[18] X. Ji, Nucl. Phys. B402 (1993) 217. 

[19] M. Maul, B. Ehrnsperger, E. Stein, A. Schafer, Z. Phys. A356 (1997) 443. 

[20] A.H. Guth and D.E. Soper, Phys.Rev. D12 (1975) 1143. 

[21] For a uniqueness theorem see, E. Carlson, Thesis, Univ. Uppsala, 1914; 

E. C. Titchmarsh, Theory of Functions, (Oxford University Press, Oxford, 1939), Chapt. 9.5. 

[22] R.L. Graham, D.E. Knuth, and O. Patashnik, Concrete Mathematics, 2nd edition, 
(Addison- Wesley, Reading/MA, 1994). 

[23] J.L. Miramontes and J. Sanchez Guillen, Z. Phys. C41 (1988) 247. 

[24] K. Bitar, PW. Johnson, and W.-K. Tung, Phys. Lett. B83 (1979) 114. 

[25] A. Dcvoto, DW. Duke, J.F. Owens, and R.G. Roberts, Phys. Rev. D27 (1983) 508. 

[26] F. Eisele, M. Gluck, E. Hoffmann, and E. Reya, Phys. Rev. D26 (1982) 41; 
R.M. Barnett, Phys. Rev. Lett. 48 (1982) 1657; Phys. Rev. D27 (1983) 98; 
R.G. Arnold et al, Phys. Rev. Lett. 52 (1984) 727; 

J.F. Gunion, P. Nason, and R. Blankenbecler, Phys. Rev. D29 (1984) 2491. 
[27] K. Varvell et al., BEBC collaboration, Z. Phys. C36 (1987) 1. 
[28] J.M. Cornwall and R.E. Norton, Phys. Rev. 177 (1969) 2584. 



33 



[29] H. Burkhardt and W.N. Cottingham, Ann. Physics (New York) 56 (1970) 453. 

[30] C.G. Callan and D.J. Gross, Phys. Rev. Lett. 22 (1969) 156. 

[31] A.V. Efremov, O.V. Teryaev and E. Leader, Phys. Rev. D55 (1997) 4307. 

[32] R.L. Jaffe and X. Ji, Nucl. Phys. B375 (1992) 527. 

[33] R.M. Barnett, Phys. Rev. D14 (1976) 70; 

R.K. Ellis, G. Parisi, R. Petronzio, Phys. Lett. B64 (1976) 97; 

R. Barbieri, J. Ellis, M.K. Gaillard, and G.G. Ross, Phys. Lett. B64 (1976) 171; Nucl. 

Phys. B117 (1976) 80; 

R. Brock, Phys. Rev. Lett. 44 (1980) 1027. 



34 



